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Abstract 

In this paper, we propose a method to obtain a constrained approximation 
of a rational Bezier curve by a polynomial Bezier curve. This problem is 
reformulated as an approximation problem between two polynomial Bezier 
curves and is solved by linear least-squares methods. The efficiency of the 
proposed method is tested using some examples. 
Keywords: Rational Bezier curves, Approximation, Bezier curves, 
Least-squares method 

1. Introduction 

Rational Bezier curves play a significant role in Computer Aided Design 
systems. However, because the forms of derivatives are quite complex and 
integrals do not exist for high-degree rational Bezier curves, the problem 
of approximating rational functions with polynomials has been raised and 
studied. Sederberg and Kakimoto pQ presented the hybrid polynomial ap- 
proximation to rational curves for the first time in 1991. Wang et al. [2] 
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presented Hermite polynomial approximations to rational Bezier curves and 
investigated the convergence condition for the polynomial approximation of 
rational functions and rational curves. Floater j3] constructed a high-order 
approximation of rational curves using polynomial curves. Huang [3] pre- 
sented a simple method for approximating a rational Bezier curve with a 
Bezier curve sequence based on degree elevation. Lee [5] converted a poly- 
nomial approximation of a rational Bezier curve to an approximation about 
two polynomial curves and obtained an approximate solution by solving a 
linear least-squares problem. Recently, the sample-based polynomial ap- 
proximation of rational Bezier curves and approximating conic sections by 
constrained Bezier curves of arbitrary degree were investigated by Lu [B] and 
Hu [7], respectively. 

The above-mentioned methods have their advantages and disadvantages. 
The conversion of a polynomial approximation of a rational Bezier curve 
to an approximation between two polynomial Bzier curves and solving the 
problem based on a least-squares method has the following advantages: First, 
the degree of approximation to the polynomial curve is arbitrary; second, for 
polynomial curves, the integral expression based on a least-squares method 
is analytic; third, because there exist a variety of methods for solving linear 
equations, computing the control points of a Bezier curve is quick and stable. 
In this paper, we extend the method presented in paper [S] to the constrained 
approximation of a rational Bezier curve. Our method also differs from that 
presented in [5] in that we don't use a degree elevation formula for Bezier 
curves. 

The paper is structured as follows. Section 2 presents some basic concepts 
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and properties regarding the problem of the constrained Bezier approxima- 
tion of a rational Bezier curve. Section 3 brings an complete solution to the 
problem formulated above in the L 2 norm. Section 4 presents some numerical 
examples to verify the accuracy and effectiveness of the method. 

2. Preliminaries 

2.1. Definitions and Properties 

A standard-form nth degree rational Bezier curve is defined as follows: 

n 

x(t) E^rw 

8=0 

where Bf(t) = — t)™ - ' are the Bernstein polynomials, p i G W 1 are 

control points and Ui G IR + (cjo = u> n = 1) are the scalar weights. 

The following are pertinent properties used in this paper: 

Property 1. The product of two Bernstein polynomials and its integral sat- 
isfies [8] 

mt)B?(t) = U^be**®, ( 2 ) 



/ B?(t)B?dt= (3) 
Jo (m + n + 1) ( ) 

Property 2. Let the two polynomials f(t) and g(t) of degree m and n with 
coefficients f™ and gf be as follows [9] 

m n 

fit) = frsm, g(t) = 5>r*rw- 

1=0 8=0 
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Their product is a degree m + n polynomial 



f(t)g(t) = £ 



m+n 



i=0 






B™ +n (t). 



(4) 



Property 3. For each appropriate function f{x), there is a unique least- 
squares polynomial approximation of degree at most n that minimizes the 
function [II Oj . 

2.2. Statement of the approximation problem 

The problem of approximating rational Bezier curves in (1) by constrained 
Bezier curves of arbitrary degree is that of finding control points q , q 1: q m , 
which define a Bezier curve of degree m 



such that the following two conditions are satisfied simultaneously [7]: 

1) Bezier curve Q(t) has C°, C 1 or C 2 continuity at both endpoints of the 
rational Bezier curve P(t). 

2) A distance function d(P, Q) between P(t) and Q(t) in the L 2 norm 
is minimized. 

In general, the integral of a rational Bezier curve over interval [0, 1] is 
either very complex or is no analytical expression because the degree of n 
grows to be large; thus, analytical solutions of Bezier curve Q(t) cannot be 
directly obtained by the least-squares method. On the other hand, although 
a rational Bezier curve cannot expressed by a polynomial precisely, it can be 
approximately represented as a polynomial Bezier curve of sufficiently high 



in 




(5) 
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degree m. Thus, we can convert the polynomial approximation of a rational 
Bezier curve to an approximation between two polynomial curves. That is, 
according to equation (1) and (5), we have 

P(0 = ^«Q(*)^*(*)«w(t)Q(«)- 

Let 

y{t)=u(t)Q(t), 

and condition 2) can be modified to the following expression 

2a) A distance function d(x, y) between x(t) and y(t) in the L 2 norm is 
minimized, that is 



i 

2 



min d(x, y) = || x(t) — y(t) \\ dt 



(6) 



x(t) -u(t)Q(t) |r dt 



3. Algorithm for constrained Bezier approximation 

The approximation algorithm is completed through a two-step process. 
First, we calculate the constrained control points of the approximation Bezier 
curve Q(t) expressed in (6). Next, we calculate the others in terms of the 
least-squares method. 

3.1. Constrained conditions of the approximation curve 
According to equation Q, we have 

x(t) = P(t)w(t). 
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Taking the derivative of the above equation r(= 0, 1, 2) times, we obtain 

3=0 



Thus, we have the zero, first and second derivatives of P(t) at two endpoints 
(t = 0, 1), satisfying the following, respectively 

P(0) = p , P(l) = Pn , 
P '(0) = nu 1 (p 1 - p ), P '(1) = nu n ^(p n - p n _ x ), 

P "(0) = n {u 2 {n - l)p 2 + 2wi(l - nu 1 )p 1 + [2ui(nui - 1) + u 2 {l - n)]p } , 
P "(1) = n{[2u n -i(nu n -i - 1) + w n _ 2 (l - n)]p n + 2w n _i(l - nujn-^p^ 
+ (n- l)w„_ip„_ 2 }. 

Then, matching the function value and derivatives up to the second order 
at both endpoints of P(t) and Q(t), the constrained control points of Q(t) 
are 

Qo = P0i Qm = Pni 

if to > 2, q 1 = — u 1 p 1 + 1 u l ) p , 

m \ m / 

/ 77, \ 

if to > 3, qr m _ x =1 w n _i p n + — £*;„_!?„_!, 

\ to / to 

if to > 4, <? 2 = — 7 -{n(n — l)u; 2 p 2 + 2nwi(m — nwi)/?, 

m(m — 1) 

+ [to(to — 1) + 2nui(nuj 1 — m) + nw 2 (l — n)]p }, 
if to > 5, q m _ 2 = — — — — {[2nu n -i(nu n -i - to) + nu; n _ 2 (l - n) + to(to - 1) 

+ 2nw„„ 1 (m - nw^ijp^i + n(n - l)w„_ 2 p„_ 2 }- 
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Accordingly, Q(t) can be rewritten as 



m— h— 1 



(7) 



i=0 



i=m—h 



where ^(z = A; + 1, A; + 2, m — h — 1) are the unknown control points and 



k = < 



m = 1, 

1 m = 2,3, & 

2 m > 4, 



m = l,2, 

1 m = 3,4, 

2 m > 5. 



3.2. Unconstrained control points of the approximation curve 

To minimize the objective function d expressed by equation ([6]), the 
derivatives of d about the points q^i = k + l,...,m — h — 1) should be 



zero; that is, = 0. Thus, we have 



Q{t)u 2 (t)B] n (t)dt= / x(t)u{t)B™(t)dt. 



(8) 



Substituting (7) into fl8J), it yields 



m— ft— 1 



£ / w 2 (*)5TW5r(*)« i d*+ £ / uj\t)BT{t)BT{t)q t dt (9) 
i=0 - 70 i=fc+l - 70 

+ J] / u\t)B?{t)B?{t) qi dt= / x(t)u(t)Br(t)dt, 
u Jo Jo 

i=m—h 



By deriving equation (9) based on properties 1 and 2, we obtain 
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2n+k I min(k,i) tm\(2n\ \ 2n+m-h-l I min(m-h-l,i) /m\ / 2n \ 

y- w y- y- w 

/ , I /2m+2n\ vv l-j1i I ^ Z^ I Z^ /2m+2n\ KV »~ 

j=0 V j=max(0,i~2n) \ i+l ) I i=k+l \j=max(k+l,i-2n) \ / 

2n+m / min(m,i) /m\ I 2n \ 

^ Z^ I Z^ /2m+2n\ Vy ^-J l li 

i=m—h \j=max(m—h,i—2n) \ £+2 ' 

_2m + 2n + l^/ g) (^) 

(Z = fc + l,A; + 2,...,m-/i- 1), 



where 

min(n,i) ^n^ / n \ 
W i = E \2n? Ui-M- 

To easily evaluate the unknown points q^i = k+1, ...,m—h—l), we rearrange 
equation (10) to obtain 

m—h—l 2n (2n\ (m\ _ n n fn\ (n\ 

(2n+2m\ VV i^i 2 m + Tl + 1 2^2^ (2n+m\ U i U iPi 
i=k+l j=0 V / i=0 j=0 Vi+j+i^ 

fc 2n f2n\ /m\ m 2n (2n\ tm\ 

-EE 7^4^- E en) 

i=0 j=0 V i+i+Z / i=m-h j=0 \ i+j+l ) 

(I = k + 1, ...,m - /i - 1). 

Finally, by property 3, we conclude that the solutions of linear system 
(11) are unique. They can be obtained by any methods introduced in [TU] . 
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4. Error estimation and implementation 

For the convenience of estimation, we use maximum distance d max to 
evaluate the approximation results. d max are the sampling-based various 
distances between the rational Bezier curve P(t) and its approximation curve 
Q{t). That is, 

d max = max \\P(t) - Q(t)\\ 

0<*<1 

Example 1. Consider a conic rational Bezier curve P(t) with control 
points p Q = (0,0), p t = (0.3, 1.5), p 2 = (1,0), and weights u = = 
0.8, U2 = 1. A 4th-degree Bezier curve Q(t) approximates P(t) with the 
contact order (1, 0) of continuity at two endpoints. The control points of Q(t) 
are (0,0), (0.12, 0.6), (0.3683, 0.9771), (0.7289, 0.5986), (1, 0). The maximum 
error distance is d max = 2.1 x 10" 3 (see Figure 1). 




Figure 1: (a) A conic P(t) with control points; (b) A 4th-degree Bezier approximation 
curve Q(t) with contact order (1,0) of continuity at two endpoints and its control points; 
(c) Error distance curve 

Example 2. Consider a cubic section P(t) with control points p = 
(0,0), = (0.2, 1.5), p 2 = (0.8, 1.5), p 3 = (1,0), and weights u — 1, u> a = 
1.2, &2 — 1-5, = 1. A 6th-degree Bezier curve Q(t) approximates P(t) 
with the contact order (2,2) of continuity at two endpoints. The control 
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points of Q(t) are 

(0, 0), (0.12, 0.9), (0.3552, 1.314), (0.5429, 1.49), (0.718, 1.035), (0.85, 1.125), (1.0, 0). 
The maximum error distance is d max = 1.99 x 10~ 2 (see Figure 2). 




(a) (b) (c) 

Figure 2: (a) A cubic Pit) with control points; (b) A 6th-degree Bezier approximation 
curve Q(t) with contact order (2, 2) of continuity at two endpoints and its control points; 
(c) Error distance curve 

Example 3. Consider a 4th-degree rational Bezier curve P(t) with con- 
trol points p = (0,0), Pl = (0.2, 1.5), p 2 = (0.4, 1.7), p 3 = (0.8, 1.5), p 4 = 
(1,0) and weights ojq = 1,U\ = 1.2, 0J2 = 1.4, u 3 = 1.2, 0J4 = 1. A 6th-degree 
Bezier curve Q(t) approximates P(t) with the contact order (2,2) of conti- 
nuity at two endpoints. The control points of Q(t) are 

(0, 0), (0.16, 1.2), (0.3008, 1.528), (0.4508, 1.5738), (0.5872, 1.528), (0.84, 1.12),(1, 0). 
The maximum error distance is d max = 1.7 x 10~ 3 (see Figure 3). 

Example 4. Consider a 4th-degree rational Bezier curve P(t) with con- 
trol points Po = (0,0), Pl = (0.2, 1.5), p 2 = (0.5, 1.0), p 3 = (0.8, 1.5), p 4 = 
(1,0) and weights uo — 1, 0J\ — 1.2, 0J2 = 1.4, uj 3 = 1.2,(^4 = 1. A 6th-degree 
Bezier curve Q(t) approximates P(t) with the contact order (1, 1) of conti- 
nuity at two endpoints. The control points of Q(t) are 

(0, 0), (0.16, 1.2), (0.3526, 1.1270), (0.5001, 1.2372), (0.6475, 1.127), (0.84, 1.2), (1, 0). 
The maximum error distance is d max = 3.4942 x 10 -4 (see Figure 4). 
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(a) (b) (c) 

Figure 3: (a) A 4th-degree rational Bezier curve P(t) with control points; (b) A 6th-degree 
Bezier approximation curve Q(t) with contact order (2, 2) of continuity at two endpoints 
and its control points; (c) Error distance curve 



(a) (b) (c) 

Figure 4: (a) A 4th-degree rational Bezier curve P(t) with control points; (b) A 6th-degree 
Bezier approximation curve Q(t) with contact order (1, 1) of continuity at two endpoints 
and its control points; (c) Error distance curve. 
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5. Conclusion and future work 

In this paper, we have proposed an algorithm for approximating rational 
Bezier curves. Compared with other methods, our method is simple, linear 
and stable. In the future, we will extend this method to solve NURBS 
approximation problems. 
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